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Abstract. The accepted explanation of the observed dichotomy of two orders of magnitude between in situ measurements of 
electron density in Hn regions, derived from emission line ratios, and average measurements based on integrated emission 
measure, is the inhomogeneity of the ionized medium. This is expressed as a "filling factor", the volume ratio of dense to 
tenuous gas, measured with values of order 1(T 3 . Implicit in the filling factor model as normally used, is the assumption that 
^ • the clumps of dense gas are optically thin to ionizing radiation. Here we explore implications of assuming the contrary: that 

the clumps are optically thick. A first consequence is the presence within H n regions of a major fraction of neutral hydrogen. 
C~-». , We estimate the mean H°/H + ratio for a population of Hii regions in the spiral galaxy NGC 1530 to be the order of 10, and 

support this inference using dynamical arguments. The optically thick clumpy models allow a significant fraction of the photons 
generated by the ionizing stars to escape from their H n region. We show, by comparing model predictions with observations, 
that these models give an account at least as good as, and probably better than that of conventional models, of the radial surface 
brightness distribution and of selected spectral line diagnostics for physical conditions within H n regions. These models explain 
how an H n region can appear, from its line ratios, to be ionization bounded, yet permit a major fraction of its ionizing photons 
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. to escape. 

^: 

q , Key words. ISM: general - ISM: HII regions - ISM: clouds - Galaxies: ISM - Methods: numerical • 

03 ■ 1. Introduction: line emission diagnostics in Hll al. 1983), but for Hii regions in Milky Way sized galaxies with 

! regions near solar metallicities either detailed photo-ionization model 

? ■ fitting, or semi-empirical methods using bright lines are used. 

^ ; Emission lines from different ionization states of the most ne of the most used of these methods involves a calibrated 
^ ! abundant elements are one of the most valuable and practical p i ot of [On]A3727A+[Om]^4959A+5007A normalized by 
tools for determining element abundances in galaxies. They are the n/3 emission line strength and plotted versus the oxygen 
relatively easy to detect and quantify. The lines emitted by H n abundance O/H, proposed by Pagel et al. ( 1 979), and subject to 
regions are intrinsically bright and sharp compared with stellar a critical review by Vflchez (1988). This plot, and similar plots 
absorption features, so easy enough to pick out and measure. for N and S , can be used in conjunction with an estimate of the 
If the electron temperature (T e ) within an Hn region can be intrinsic radiation hardness from the ionizing stars: + ,0 ++ , 
determined, it is not too difficult to determine the overall abun- S +/ S ++ (Vflchez & Pagel 1988), to measure abundances and 
dance of a given element by measuring the relative emission the abundance gradient across the face of a spiral galaxy. An in- 
line strengths of lines due to several ions. To estimate T e , the teresting early summary of oxygen abundance gradients in spi- 
excited [Om]/i4363A line can be measured and its ratio with ra i s obtained using these principles can be found in Edmunds 
the /14959A and 5007 A transitions to the ground state, yields (1989). 
a good r e value. This works well at intermediate metallicities 

(e.g. in the Magellanic clouds) but not at near solar metallic- As early as 1967 Peimbert pointed out that temperature 

ities, where the lower Hn region temperatures lead to a weak fluctuations in Hn regions will cause temperature estimates 

lambda 4363A line, or at very low metallicities where all the ma de on the basis of forbidden line ratios to differ from those 

Om lines are weak. For local Hn regions T e can be deter- made via the ratio of Balmer lines to Balmer continuum, and as 

mined well from radio recombination lines (see e.g. Shaver et sucn will lead to imprecision in any physical parameter deter- 

mination requiring electron temperatures, notably that of abun- 

Send offprint requests to: C. Giammanco dances. Rubin (1969), considered explicitly the effects of tem- 
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perature fluctuations on abundance determinations, showing 
that for certain excitation conditions, neglect of these fluctua- 
tions will lead to underestimates of the abundances of elements 
such as oxygen and nitrogen. Rubin (1968), also explored the 
effects of density fluctuations, finding that the smaller and 
denser are the clumps of dense gas within a region, the higher 
the resulting effective electron temperature. He calculated mod- 
els whose volume emissivities in Ha and [Om]/15007A are af- 
fected by density and its fluctuation, but these were of the kind 
represented in the "filling factor" approximation which we will 
consider in some detail below, and as such are not an especially 
realistic representation of the density distribution within an H n 
region. In practice, the effect of density fluctuations on the ion- 
ization equilibrium of a gaseous nebula has been considered 
essentially only as the cause of temperature fluctuations. This 
was clear in early work on the Orion Nebula by Osterbrock & 
Flather (1959) and in a subsequent more detailed examination 
of the same H n region by Simpson (1973), was brought out in 
abundance determinations in the Magellanic Clouds by Dufour 
(1975), and also noticed by Krabbe & Copetti (2002) when 
considering the effect of the high density knots in 30 Doradus 
in slightly reducing the local electron temperature in apparent 
contradiction with Rubin (1968). In all of these studies, the un- 
derlying assumption is that there is a continuously varying am- 
plitude spectrum of fluctuation size and density, although since 
the emission intensity of a recombination line scales as the 
square of the electron density, any diagnostic diagram will give 
the strongest weighting to the densest fluctuations. However, 
relatively little direct attention has been given to the effects 
of the density fluctuations themselves on the radiative trans- 
fer properties of an Hn region. Roy et al. (1989) attributed the 
variation of a number of line ratios on small scales within a 
number of large luminous extragalactic H n regions, to small 
scale density variations within each region, but no direct atten- 
tion has been paid to the possible effects of major density fluc- 
tuations on the overall energy budget of a region. In the present 
paper we address this problem, using a straightforward mod- 
elling technique to represent the presence of dumpiness within 
the regions. The clumps in the model show considerably more 
contrast with their surroundings than is assumed for a continu- 
ous density variation distribution, but the model has a density 
distribution much more like that empirically observed in the lo- 
cal neighbourhood (the only place where individual ~ 1 pc scale 
clumps in the ISM can be separately detected and measured). 
These clumps are also in better physical accord with the two 
and three phase models of the ISM which have been proposed 
to explain its observed properties within the Galaxy (Cox & 
Smith 1974; McKee & Ostriker 1977). 

In the following section we give an explanation of the 
model and its assumptions, and in subsequent sections we dis- 
cuss the basic structure of the ionization distribution within the 
model. In Sect. 4 we compare the predictions of the generic 
radial surface brightness distribution of an individual region 
with the observationally derived distributions for Hn regions 
in external galaxies. Comparison of some predicted emission 
line ratios for this and previous models with observed line ra- 
tios, are discussed in Sect. 5, and finally, in Sect. 6 we analyze 
the important implications of this model for the observationally 



derived masses of gas within Hn regions. In the conclusions 
we bring out in what ways these clumpy models are qualita- 
tively different from both the classical Stromgren model with 
homogeneous gas, and its derivatives obtained using the stan- 
dard "filling factor" approach. 

2. Inhomogeneities: the classical filling factor or 
optically thick clumps? 

The possible effects of density inhomogeneities on the propa- 
gation of ionizing radiation in H n regions were first discussed 
in some detail by Stromgren (1948). However the standard 
way to take into account the fact that the ISM is highly non- 
homogeneous on scales of parsecs is to scale the optical paths 
and the volume emissivity by a filling factor, following the ini- 
tial prescription of Osterbrock & Flather (1959), which we will 
term here the FF approximation. The technique assigns a frac- 
tion of the total volume of an H n region to relatively dense gas, 
with the remaining volume taken to be of negligible density. 
The fraction occupied by dense gas is termed the filling factor, 
and it takes values of order 10~ 3 in typical observed cases. The 
need to treat the volume of a region in this way is because the 
measurement of the in situ electron density using line ratios ei- 
ther at radio wavelengths (e.g. Dravskikh & Dravskikh 1990) 
or optical wavelengths (e.g. Copetti et al. 2000) gives values of 
order a few times 10 2 cm 4 , whereas the mean electron density 
obtained from the emission measure in Ha is of order 1-2 ctrr 3 
(Rozas et al. 1996). In the form introduced by Osterbrock & 
Flather (1959) and in general use for dealing with the transfer 
of Lyman photons in H n regions, the underlying assumption is 
that the dense clumps are small enough to be optically thin, so 
that each clump is fully ionized. The overall effect is a statis- 
tical modification of the classical Stromgren sphere structure, 
in which radiation from the central ionizing stars ionizes mate- 
rial within a given radial distance from the stellar cluster and is 
effectively fully absorbed within the H n region, except for the 
minority of cases where either the placental cloud is sufficiently 
small for some radiation to escape, or the cluster is formed ec- 
centrically within the cloud, and thus near the cloud edge in one 
direction. Although the physical conditions for recombination 
line formation are biased towards the high densities, since these 
emission line strengths are proportional to the square of the 
electron density, and hence this model in which only the dense 
clumps are considered, gives many valid results, we believe that 
a better description of the transfer is given by using a medium 
with clumps whose optical depth is high. Rather interesting ev- 
idence for this is found within the local ISM where a detailed 
examination of structure is possible within a zone which lies 
within a major H n region based on the OB associations towards 
Scorpio and Centaurus. Here Trapero et al. (1992, 1993) found 
that more than half the mass of the ISM within 300 pc of the 
Sun is in the form of dense (~ 100 ctrT 3 ) compact clouds with 
characteristic sizes of order a few parsecs, while the remain- 
der forms a warmer substrate with densities lower than 1 crrT 3 . 
Inspection of optical images of H n regions within the Galaxy, 
e.g. the Orion Nebula, confirms this picture of ionized mate- 
rial in dense blobs or clumps several parsecs across, separated 
by larger volumes of less dense gas. In fact, this structure con- 
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forms to the rather clean separation into cool and warm phases 
first proposed by Field, Goldsmith & Habing (1969) and later 
extended to three phase models which include a hot component 
by Cox & Smith (1974) and by McKee & Ostriker (1977). A 
particularly clear theoretical underpinning of these multi phase 
models is given in the classical text by Spitzer (1978), who 
explains in terms of detailed balance, why certain ranges of 
density and temperature in the ISM are unstable, and how this 
leads to a geometrical differentiation of the stable phases. The 
transfer properties of an H n region with a clumpy structure in 
which there is a clear phase separation, with correspondingly 
large well defined clumps are significantly different from those 
of a traditionally modelled clumpy region where the structure is 
not scaled, and it is this difference which we explore here. We 
will term the traditional model the FF model, and our model 
the "clumpy" model, since the density and volume contrasts in 
our model are greater. The clumps in our model, and in the H n 
regions themselves, are optically thick to the Lyman contin- 
uum (Lyc) and this has significant effects, e.g. implies that the 
Stromgren sphere approximation loses validity, and that tradi- 
tional diagnostics for photon escape such as those applied in 
McCall, Rybski & Shields (1985) must be re-quantified. 

To give an analytical first order approach which is easy to 
follow we first take our clumps to be uniform spheres, with a 
characteristic size d and a cross-section for radiation intercep- 
tion cr. A clump situated at a distance R from the centre of a 
(spherically symmetric) region intercepts a fraction crj'XnR 2 of 
the ionizing radiation, and for a number density of clumps n, 
the clumps in a shell of thickness 6 between R and R + 6 from 
the centre, intercept a fraction of the radiation AnR 2 5n(cr / '4nR 2 ) 
i.e. drier. We term this fraction /, which represents the proba- 
bility that a photon is absorbed within this shell, i.e. the fraction 
of photons absorbed within the shell. With the simplifying as- 
sumption that the number density of clumps is not so large that 
there is a significant probability that a ray from the centre of 
the region cuts more than one clump on its path to the surface, 
the probability e for a photon to escape from the region is given 
by 

e = d-ff (1) 

where N is the total number of shells. We are particularly inter- 
ested in the production of Ha by recombination in a process 
which can be considered as a down-conversion of Lyc pho- 
tons, terming the intensity of Ha emission by the f th shell 
from the centre as /,. We can formulate this in terms of the 
intensity, If, which would be produced within a filled shell 
which has a thickness d, multiplied by the fraction of the shell 
which absorbs photons, i.e. the fractional area occupied by 
dense clumps, giving 

I i = l{ a f(l-ff-» (2) 

Both equations, ^and|2] use a thin spherical shell approxi- 
mation, which will not be valid near the centre of the region, but 
in any case, if n is small the contribution of that shell to the to- 
tal emission can be neglected, not only for simple geometrical 
reasons, but because in the innermost few parsecs the clumps 
will be dissipated by a combination of stellar winds and radi- 
ation. For simplicity we will set the shell thickness as the size 
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of a single clump, and we can thus neglect the innermost shell 
in any calculations. The fractional error in doing this for a shell 
where the clumps were not dissipated would be just /, (and 
we will see that this ranges from 0.001 to 0.01 so this would 
be no problem) but in practice the error will be smaller than 
/. Even if our computation were incorrect for the first three or 
four shells, the error implied would be small. In particular, if 
the clumps in the second shell were either dissipated or fully 
ionized, our error would again be of order /, and this problem 
will not be important, even in regions with the most luminous 
ionizing stars, for shells beyond the third or at most the fourth, 
for clumps of scale size ~ 1 pc. 

It is instructive to use the spherical clump approximation 
to derive some quantitative relations associated with the filling 
factor, i.e. the fractional volume occupied by dense gas. For 
spherical clumps of radius r, for which cr is nr 2 we find 

/ = n?n (3) 

and setting the shell thickness delta to be the clump diameter 
(6 = 2r) we obtain 

/ = (3/2)v c » (4) 

where v c is the volume of a clump. Let us define a geometrical 
filling factor <pc as the fraction of the total volume of the H n 
region occupied by the clumps and so we find 

<Pc = (4/3)r//<5 (5) 

and if we put 5 - 2r we have 

/ = (3/2)00 (6) 

It is important to be explicit here about the method nor- 
mally used to infer filling factors of dense gas in Hn regions. 
Following Osterbrock (1989, p. 153), the luminosity L emitted 
by an H n region whose dense clumps have a total volume V c 
and electron density A^ e , in a given spectral line is 

L = kN e 2 V c (7) 

where A; is a constant which depends only on which emission 
line is being observed, and not (to first order) on the physical 
conditions. If we assume that we are dealing with an Hn region 
whose total volume is V, then the ratio V c /V is the filling factor, 
<p so we can express Eq.0as 

L = kN e 2 4> V (8) 

which for the case of Ha enables us in principle to determine 
the filling factor <\> Q using 

O = k'LHjdAnl^N 2 ) (9) 

where R is the radius of the Hn region, and the constant k' 
takes the value appropriate for the Ha emission line. We note 
especially that in an FF model the two filling factors, <pc and <p a 
are equal, since all clumps are fully ionized, but in the present 
clumpy model this is not the case. We will see that in general 
<p < §g since most clumps are not fully ionized, and the differ- 
ence between these quantities can be easily an order of magni- 
tude. In addition, deriving <f> using an integrated emission line 
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Fig. 1. Equatorial section of a spherical clump with density 100 
cm -3 and radius 1 pc illuminated by an ionizing source and 
placed at difFerent distances from it. The light area represents 
the ionized volume of the clump (seen in cross-section) and the 
dark area represents the volume which remains neutral. Gray 
areas show partial ionization. In this example, obtained with 
CLOUDY, we used an ionizing source emitting 10 42 erg s _1 
yielding conservatively high clump ionization. 



luminosity and Eq. ^implies obtaining a reliable value for the 
radius of the H n region, which is not the most accurate of ob- 
servations because the projected surface brightness of a region 
blends into the diffuse background, so that any measurement 
which relies on establishing a clear boundary entails errors. We 
will discuss these points further, later in the paper. 

3. Calculations using FF and clumpy models 

In order to illustrate both qualitatively and quantitatively the 
differences which arise using the two types of models we have 
generated specific models to compute line strengths and ra- 
tios under defined conditions. All these models rely on the 
CLOUDY suite of programs (see Ferland et al. 1998). For the 
FF models we have used CLOUDY in a standard way. We have 
generated a tree family of models of different metallicity with 
a set of stars: 3, 10, and 30 03 stars at the centre, a mean gas 
density of 100 cm 4 and filling factors of 10~ 3 , 10~ 4 , values 
designed to match the optical filling factors obtained observa- 
tionally (e.g. Rozas et al. 1996). The outer radius of the region 
is a parameter which can be varied, so that if we start with 
low values and progress to high values we go from a density 
bounded to an ionization bounded regime. We used essentially 
solar metallicities, as specified in the manual of CLOUDY, and 
for the purpose of streamlining the calculation we simulated an 
03 star by a black body with a temperature of 5 1230 K, and an 
ionizing flux 10 39 31 erg s _1 (Vacca et al. 1996). 

Since the clumpy models are based on clumps of denser 
gas whose properties: size and gas number density, are defined 
as inputs, we first made a set of simple tests in which a clump 
of radius 1 pc and density 100 cur 3 was illuminated by an 
ionizing source at different distances. The source was given an 
ionizing luminosity of 10 52 5 photons s _1 which corresponds to 
almost 300 03. This is a very luminous source, so that our con- 
clusions about the ionized fractions of the clumps at varying 
distances from the source are conservative. We can see the re- 
sult of this test in Fig.^ where the fractional volume ionized, 
seen in cross-section, and the degree of ionization, are shown 
as functions of the distance of the clump from the source. 



It is obvious at a glance that at distances further from the 
source than some 10 pc, a clump is essentially opaque to the 
ionizing radiation, and in terms of simple transfer can be vir- 
tually considered a black disc. However with the filling factors 
found observationally, the net effect of a clumpy region where 
the gas is concentrated into these clumps will be to act as a "ra- 
diation sieve" allowing a major fraction of the ionizing photons 
to pass between the clumps and escape. It is also clear that at 
distances beyond some 20 pc from the source the major frac- 
tion of a clump is shielded and left unionized. Depending on 
the details of the distribution of the clumps within the H n re- 
gion, the net effect of integrating unionized and ionized clump 
masses will be that a relatively small fraction of all the clump 
mass in the H n region will be ionized. It is this fraction which 
is detected as the "optical" filling factor. A useful canonical es- 
timate of the ratio of neutral to ionized gas in an H n region of 
radius 150 pc, i.e. one of the ultraluminous "giant" H n regions, 
is a factor 10, as we will show empirically in section 6 below, so 
that in our clumpy models we use a "geometrical" filling factor 
of 10~ 2 , 10~ 3 (including both H° and H + ) which in broad terms 
corresponds to an optical filling factor of 10 3 , 10 ~ 4 (including 
H + only). 

In the computations made using the clumpy models we 
have used the formalism implied in Eq.[2]in conjunction with 
the CLOUDY program of continuum and line transfer (Ferland 
et al. 1998). 

4. Differences in the radial distributions of ionizing 
and recombination photons between the FF and 
clumpy models 

A key qualitative difference between FF models and clumpy 
models is the prediction of the fraction of ionizing photons 
which escape from an Hn region. In an FF model, as in a 
classical Stromgren uniform model, there is a clear distinction 
between density bounding and ionization bounding, which de- 
pends only on the mean gas density and the modelled radius of 
the region. In slightly simplified terms once a region reaches 
a specified size for its density distribution it will be ioniza- 
tion bounded, so no photons escape, while for smaller radii 
it will be density bounded. In a clumpy model this clear dis- 
tinction is not maintained, since a fraction of ionizing photons 
escapes from essentially all clumpy models with reasonable 
physical parameters, and the only question is to quantify that 
fraction (Eq. [Q. In other words, all regions modelled using a 
clumpy model will be, de facto, density bounded. However, in 
each individual clump (except for the minority of clumps close 
to the ionizing source) all photons incident on the clump are 
absorbed. Thus the emission line ratios for lines produced by 
down-conversion within the clump will be those for ionization 
bounded regions, and as essentially all the line emission from 
within a region comes from the clumps, any diagnostic diagram 
using line ratios will give an ionization bounded signature. In 
other words clumpy H n regions will always be diagnosed as 
ionization bounded, although in terms of photon escape they 
behave as density bounded. An illustration of these different 
types of behaviour is given in Fig.[2]where we have plotted es- 
cape fraction against H n region radius for three FF models and 
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Fig. 2. Dependence of photon escape fraction on the H n region 
radius for different models. In classical FF models, the escape 
fraction depends on the radius, on the filling factor, and on the 
strength of the ionizing radiation (we have considered ionizing 
clusters with 3, 10, and 30 03 type stars, and a filling factor 
of 10~ 4 ). In contrast, in the clumpy models the dependence is 
only geometrical; in other words, the escape fraction depends 
on the radius of the region, and on the radius and density of the 
clumps. We show results of three clumpy models (labeled with 
C) assuming different filling factors 10~ 2 , 10~ 3 and clump radii 
(1 and 2 pc). 



Fig. 3. Observationally derived Ha radial emissivity profile of 
a bright H n region of the barred spiral galaxy NGC 1530 (solid 
line). The profile has been normalized to facilitate its compari- 
son with the profiles predicted by clumpy and FF models. The 
spatial resolution of of the four theoretical models has been de- 
graded to match the resolution of the Ha image of NGC 1530 
(approx. 50 pc per pixel). We can see that the clumpy mod- 
els produce concave profiles, similar to the observations. On 
the contrary, FF models, give rise to convex profiles. The dot- 
ted line represents an ideal Stromgren model, whose ionization 
fraction is constant and equal to 1 inside the H n region volume, 
except close to the boundary of the Stromgren sphere. 



three clumpy models. The difference is clear. The FF models 
show behaviour akin to that of a classical Stromgren sphere. 
The fraction of photons escaping falls quite sharply as the ra- 
dius of the H n region reaches a critical value, which depends 
on the luminosity of the ionizing stars and the mean density 
of the gas (in the models illustrated here this latter parameter 
was maintained constant so the effect is not seen). The clumpy 
models, labeled with C, show quite different behaviour, which 
is essentially geometrically dominated. The escape fraction for 
an H ii region for a given source luminosity falls off steadily 
as the radius of the region rises, but there is no cut-off radius, 
as the fall-off is determined basically by the geometrical cover 
factor of the clumps, each of which is optically thick. In fact the 
plots for these models slightly underestimate the escape factor, 
as they have assumed no overlapping between clump cross sec- 
tions for photon blocking. 

As a direct observational test of the differences in the pre- 
dictions of FF and clumpy models we have looked at radial 
profiles in surface brightness of H n regions (Rozas et al. 1998), 
a test which makes best use of the spherical symmetry incor- 
porated into our simplest clumpy models. We selected one of 
the largest isolated H n regions with large radii in the galaxy 
NGC 1530 for which we had continuum subtracted Ha im- 
ages with good signal to noise ratio, measured its radial sur- 
face brightness profile in Ha, and from this profile derived the 
distribution of volume emissivity in Ha as a function of radius 
within the region. The radial plot of this contribution function 
is shown in Fig. [3] In the same figure we give the predictions of 
four models, for which the luminosity of the central source was 



scaled to yield the observed Ha luminosity, to a first approxi- 
mation: two FF models with different assumed filling factors 
and two clumpy models also with different geometrical filling 
factors. It is clear from the figure that the FF models give con- 
tribution functions which do not fit those of the observed re- 
gions. Although the radial profiles of these models differ de- 
pending on the filling factor, they have in common a convex 
form, with a rather sharp fall off towards the edge of the region. 
The clumpy models, on the other hand, show convex profiles, 
and give far better agreement with those derived from the ob- 
servations. The Ha volume emissivity for a uniform Stromgren 
sphere is included for comparison. 

5. Predictions of line strengths and ratios: 
diagnostic diagrams for density bounding 

A classical test for density bounding in H n regions is the plot 
of [Ora]/[Oii] vs. ([Oiii] + [Oii])/Hj0, proposed by McCall et al. 
(1985). Examination of this plot for a set of extragalactic Hn 
regions, led them conclude that most H n regions are ionization 
bounded. Other authors (Roy & Walsh 1997; Kennicutt et al. 
2000) have reached the same conclusion based on similar anal- 
ysis. However, there is good evidence (Ferguson et al. 1996; 
Zurita et al. 2000; Relano et al. 2002) that a major fraction, of 
order half, of the ionizing photons are escaping from the H n re- 
gions of external galaxies, and ionizing their diffuse interstellar 
media. This evidence is based on two types of measurements. 
Firstly the integrated diffuse Ha measured over the full face of 
a galaxy is found to depend linearly on the integrated Ha lu- 
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Fig. 4. Density bounding diagnostic diagram (McCall et al. 1985). Solid lines delimit the area occupied by McCall's data and 
the dotted line indicates where this limit is least certain. Data points correspond to H n region observations, which have been 
obtained from Cedres & Cepa (2002), Van Zee et al. (1998), Kennicutt & Garnett (1996), Walsh & Roy (1997) and Roy & Walsh 
(1997). Colour symbols refer to models (FF in a. and clumpy in b.) assuming 3, 10 and 30 ionizing 03 type stars and 3 different 
abundances: 0.1 solar, 0.5 solar, and solar. Colours indicate the escape fraction of ionizing photons in each case. We can see that 
all models, even those with a strong escape fraction, are located in the areas covered by the observational data. 
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Fig. 5. Density bounding diagnostic diagram: [Oi]/16300A vs. [On]/13727A (cf. Iglesias-Paramo & Munoz-Tunon 2002). a.) 
Sequence of FF models obtained for an optical filling factor equal to 10~ 4 . b.) Clumpy models for / = 10~ 3 . c.) Same as a. but 
for <p = 10~ 3 . d.) Same as b. but for / = 10~ 2 . The plots show that as the radius increases, the data points move towards the 
top-right region of the plots (the "ionization-bounded" area) and the escape fraction decreases. For clumpy models, Hn regions 
with high escape fraction (see Fig |5j) for / = 10~ 3 ) occupy the same general parameter space as models with moderate escape 
fraction (Fig[3J> for for / = 10~ 2 ). For the convertion of / to the geometrical filling factor <p c see Eqs.|5]and|6](/ is of order <p c ). 
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minosity of the most luminous H n regions (above a specified 
lower limit in Ha luminosity) in that galaxy (Beckman & Zurita 
2004), suggesting very strongly that escaping Lyman contin- 
uum photons from the H n regions are responsible for ionizing 
the diffuse medium, and by inference are escaping from the 
H ii regions themselves. Secondly a careful model of a specific 
galaxy (NGC 157) in which all the Hn regions, with measured 
Ha luminosities and positions were taken as sources of ionizing 
photons, which were allowed to propagate across the galaxy, 
gave excellent agreement between the modelled geometrical 
distribution of the diffuse Her emission from the whole galaxy 
and the observed distribution (Zurita et al. 2002). The fraction 
of photons which need to escape is at least 30% of those emit- 
ted by the OB stars within the regions. We have predicted some 
of the line ratios previously used for density bounding diagnos- 
tics using our clumpy models and compared these predictions 
with those of the classical FF models. In Fig. 0] we show the 
large data set presented by several authors (including McCall 
et al. 1985) in the log([Oiii]/[Oii])-log(([Om]+[Oii])/HyS) pa- 
rameter space, together with the predictions of the FF ©) and 
clumpy l|4j>) models for different metallicities. We can see that 
in both cases, the observations can be reproduced by mod- 
els with strong ionizing photon leakage (indicated by color 
scale) and models with negligible escape of ionizing photons. 
We can also infer two types of conclusions from this compari- 
son. Firstly, when we inspect in detail the comparisons between 
models and observations, the clumpy models do give a signifi- 
cantly improved account of the observations, and secondly the 
use of [O m]/[0 n] as a test for density bounding is not valid in 
either type of models. The model grids, notably for the clumpy 
models, overlap with, and reproduce the global distribution of 
the measurements. It is still true, as predicted from previous 
types of models, that a higher [O m]/[0 n] ratio tends to indicate 
a higher escape fraction of ionizing photons from a region, but 
the degree of ambiguity about this prediction for a given exam- 
ple is very high, given the error bars on a given measurement, 
for which the reader is referred to McCall et al. (1985). This 
last conclusion is in good agreement with Iglesias-Paramo & 
Munoz-Tundn (2002), who proposed an alternative diagnostic 
based on the line strengths of [Oi]/16300A and [On]A3727A. 

In Fig. |5] we compare the different predictions of clumpy 
$5$i and|5J;) and FF (|5j) and|5jl) models for the [Oi] diagnos- 
tic diagram. Different values of absorption fraction per shell 
(/ =10~ 2 and the extreme value 10~ 3 ) and optical filling fac- 
tor (0 O =1O -3 , and the extreme value 10 ) have been used in 
the clumpy and FF models respectively. Figs.|5^ and c clearly 
show that FF models with low escape fraction (dark points) 
occupy the same area of the plot independently of the filling 
factor adopted. This means that if H n regions were well de- 
scribed by FF models, this representation would be valid for 
density bounding diagnostics. However, clumpy models with 
high escape fraction in Fig. [5j>, occupy the same parameter 
space as the clumpy models with low escape fraction in Fig.|5jl. 
Therefore, before using these diagrams, we need to make as- 
sumptions about the parameter /. In order to understand this 
behaviour, we show in Fig. [6] the predicted [O i]6300A/H/3 line 
ratio as a function of H n region radius, again comparing the 
clumpy models (J6j)) with "standard" FF models ©). It is a 



general feature of all the models that as the radius of the Hn 
region increases the ratio rises, but for two different reasons: in 
the FF case, there is a rapid rise to an asymptotic value, because 
as the ionizing photons are increasingly absorbed within the re- 
gion, the fraction of hydrogen which remains in neutral state 
rises, and therefore the [Oi] emission, induced by a charge- 
exchange reaction. However the clumpy models tell a different 
story. Here there is no clear transition in the line intensity ra- 
tio because there is no clean transition in the physical proper- 
ties of the regions; in other words, each clump in which emis- 
sion lines are formed and emitted is ionization bounded itself. 
The distance of the clumps from the ionizing sources deter- 
mine different ionizing parameters, i.e. different values of the 
[O i]6300A/Hj6 line ratio which average in a complex way giv- 
ing rise to the behaviour shown in Fig.|6j). We conclude here 
that the distribution of data in McCall et al's plot is compati- 
ble with density bounding, as density and ionization bounded 
models, for both clumpy and FF, are in the same region of the 
plot. We have some reason to prefer clumps as for dynami- 
cal reasons (see next section), there does appear to be a high 
proportion of Hi within large Hn regions. However to decide 
definitively which is the best way to model Hn regions, we 
need to take into account other diagnostic diagrams which will 
be the topic of a forthcoming paper (Giammanco et al. 2004). 

6. Gas masses in H n regions 

It is very difficult to estimate the amount of neutral hydrogen 
within an H n region, because of the need to separate the contri- 
butions of hydrogen clouds along the line of sight to the object 
(and even beyond it) from the 21 cm emission from the specific 
volume of hydrogen we wish to identify. This difficulty is seri- 
ous even using Doppler techniques to try to isolate the emission 
from the region under study. For H n regions in external galax- 
ies the velocity does give us a better discriminator, but here 
angular resolution limits our ability to define the H i content 
(and also the possible H2 content) of even a large luminous H 11 
region. However there have been a few attempts to supplement 
direct measurements by indirect estimates of the neutral gas 
mass within an Hn region using dynamical arguments. Yang 
et al. (1996) estimated that in order for the basic emission line 
width (the width unaffected by locally expanding material, and 
which required careful two dimensional spectroscopy to pick 
out) of Ha from the region NGC 604 in M33 to be consistent 
with a system in virial equilibrium, the neutral gas mass had 
to be an order of magnitude higher than their observed mass 
of ionized gas. They supported their argument using an obser- 
vational HI mass estimate by Deul & van der Hulst (1987), 
and an observational estimate of the H2 mass by Viallefond 
et al. (1992), and they themselves showed that the total stellar 
mass is not a major fraction of the total mass of the region. A 
similar though less clear cut argument was produced by Chu 
& Kennicutt (1994) for 30 Doradus. Comparing emission line 
widths from many features across this huge, nearby, and lumi- 
nous H 11 region with mass estimates based on the virial theo- 
rem, they found that even an order of magnitude more H 1 than 
H n would not yield virial equilibrium, relating the discrepancy 
to dynamical energy inputs from the hot stars and supernova 
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Fig. 6. a.) Behaviour of line ratio [Oi]6300/Hj8 versus the Hn region radius for different families of Hn region models in the 
FF approximation, for two different stellar contents and metallicities. As the radius increases, the regions change from density 
to ionization bounded. The [Oi]6300/HyS ratio becomes constant when the models become ionization bounded, b.) Same as a.) 
but for clumpy models. In this case the line ratio increases very slowly because the emission lines result from the sum of the 
contribution of all the clumps, which are ionization bounded. The reason why we see a slow increase rather than a plateau is due 
to the fact that each clump has a different ionization parameter, depending on its distance from the ionizing stellar source. 



remnants. However in this case there are no corresponding lit- 
erature estimates of the neutral gas (atomic or molecular ) with 
which to compare. It is clear from the above articles that if 
we assumed that the observational estimates of the ionized gas 
gave us the complete gas mass, i.e. that there is little or no 
neutral mass within the H n regions, we would put the regions 
way out of equilibrium. It is less clear, however that they are in 
fact in equilibrium generally, though locally perturbed by stel- 
lar energy inputs. This however, was the basic result of the very 
detailed study by Yang et al. (1996) for NGC 604. 

In this context we see that we could use our clumpy mod- 
els to estimate neutral gas masses in H n regions purely on the 
basis of ionization equilibrium, once the basic postulates about 
the clump properties had been laid down. The calculation is 
straightforward. Our model allows us to estimate the fraction 
of a clump which is ionized at any chosen distance from the 
ionizing source, so once we know, or assume, the distribution 
of clumps within the region as a whole we can compute the 
neutral gas fraction. The results of calculations of this type for 
a set of observed Hn regions in the barred spiral NGC 1530 
are shown in Fig. Q The results shown are based on the spe- 
cific assumption that a constant fraction of 30% of the ionizing 
radiation from the stellar sources escapes from all the H n re- 
gions, but the qualitative result would not change if we used a 
different basis for the computation (for example increasing this 
escape factor to 50% in fact decreases the neutral/ionized gas 
mass ratio less then factor 2). The graph shows the ratio of true 
clump filling factor, </>g, which is the fractional volume of an 
H ii region occupied by dense gas clumps, to the optical filling 
factor (f> , which is the fraction occupied by the ionized portions 
of the gas clumps. We can see from Fig.0that ratios of an order 
of magnitude between the neutral gas and ionized gas masses 
are indicated in these models, as a natural consequence of the 
clumpy nature of the medium. This gives a firmer basis for the 
use of this factor in making the calculation of the gas mass in 



NGC 604, and implies that a virial equilibrium solution for the 
line widths observed by Yang et al. makes good physical sense. 
In a forthcoming paper (Relano et al. 2004, A&A, submitted) 
we will present a general study of the topic of the internal kine- 
matics of H ii regions based on observations of complete popu- 
lations of regions in the discs of spiral galaxies, in which one of 
our conclusions is that the widths, cr of the integrated Ha emis- 
sion from the majority of the H n regions observed in external 
galaxies are broader than those expected for systems in virial 
equilibrium, but that the lower limiting envelope in sigma for 
the whole population, plotted in a diagram of cr v. Ha luminos- 
ity, represents a locus of regions which are virialized.Were it 
not for the major contribution of the neutral gas in the clumps 
the regions would all appear to be gravitationally unbound. 

We can thus conclude, at least tentatively, that measure- 
ments of the internal dynamics of Hn regions, supplemented 
by direct estimates of H i in the few cases where this is possi- 
ble, give good support to the use of our clumpy models to in- 
terpret the emission parameters. This support is quantifiable in 
terms of the fractional estimates of neutral gas mass obtained 
dynamically compared with the fractional estimates based on 
clumpy models with reasonable radiative parameters. Lastly it 
is relevant to point out that Relano et al. (2004) have shown 
that the stellar mass cannot make up more than a few percent 
of the total H ii region mass for these luminous regions, and so 
does not materially affect the dynamical balance of a region. 

7. Conclusions 

The aim of this paper has been to offer an initial alterna- 
tive version to the conventional inhomogeneous models for 
H ii regions. The classical filling factor approximation assumes 
that the denser clumps are optically thin, and this gives rise 
to a statistically modified Stromgren sphere with a well de- 
fined boundary in ionization bounded conditions. The use 
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Fig. 7. Ratio between the geometrical and optical filling factors 
as a function of Ha luminosity for H n regions from the cata- 
logue of NGC 1530 (Relano et al. 2004). C has been obtained 
assuming spherical clumps with radius 1 pc and a photon es- 
cape fraction of 30%. Under this hypothesis, the mean ratio of 
filling factors is ~ 10, while assuming a photon escape fraction 
of 50% the mean ratio is ~ 6. 

of an optically thick approximation gives some substantially 
different predictions. We have shown in particular, the dif- 
ferences implied in the diagnostic diagram of [Oiii]/[On] v. 
([Om]+[On])/Hy8, a diagram used, inter alia, for separating re- 
gions which are ionization bounded from those which are den- 
sity bounded. The optically thick clumpy model gives predic- 
tions in rather better agreement with the observed data than 
the traditional filling factor models, but more critically, we 
show that in neither model is the separation of the two condi- 
tions (ionization and density bounding) very clean, and in the 
"clumpy" model there is complete overlap, rendering the dia- 
gram ineffective as a diagnostic. We apply the same exercise 
to the more recently suggested [O i]/H/? v. [O ii]/H/? diagram, 
showing that in this case the diagnostic would function ade- 
quately for the traditional models, but for the optically thick 
clumpy models the results would be ambiguous. The models 
described here, though still geometrically quite schematic, do 
appear to give an improved account of the radial profiles of 
Hn regions in surface brightness. They also predict as a nat- 
ural consequence that Hn regions should contain a majority 
fraction of H i, and this prediction seems to be borne out by 
observations of the internal dynamics of H n regions. 
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